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ABSTRACT 



Aims. We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within envi- 
ronments more comparable to loose groups, including the Milky Way's Local Group. 

Methods. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the 
field and in loose groups environments with similar properties to the Local Group. Apart from the differing environments of the 
galaxies, the samples are kept as homogeneous as possible with equivalent ranges in last major merger time, halo mass and halo spin. 
Comparison of these two samples allow for systematic differences in the simulations to be identified. A kinematic decomposition is 
employed to objectively quantify the spheroid-to-disk ratio and to isolate the disk-star population. Metallicity gradients, disk scale 
lengths, colours, magnitudes and age-velocity dispersion relations are studied for each galaxy in the suite and the strength of the link 
between these and environment of the galaxies is studied. 

Results. Metallicity gradients are consistent with observations of HII regions in spiral galaxies and, in agreement with observations, 
correlate with total galaxy mass. The bulge-to-disk ratio of the galaxies show that these galaxies are less spheroid dominated than 
many other simulated galaxies in literature with the majority of both samples being disk dominated. We find that secular evolution 
and mergers dominate the spread of morphologies and metallicity gradients with no visible differences between the two environment 
samples. In contrast with this consistency in the two samples there is tentative evidence for a systematic difference in the velocity 
dispersion-age relations of galaxies in the different environments. Loose group galaxies appear to have more discrete steps in their 
velocity dispersion-age relations, if this is true it suggests that impulsive heating is more efficient in the stars of galaxies in denser en- 
vironment than in the field. We conclude that at the current resolution of cosmological galaxy simulations field environment galaxies 
are sufficiently similar to those in loose groups to be acceptable proxies for comparison with the Milky Way provided that a similar 
assembly history is considered. 

Key words. Galaxies: Local Group, formation, evolution - methods: numerical 
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1. Introduction 

It is well established that galaxy interactions and merg- 
ers result in s ignificant ch a nges i n a system's s ta r for- 



mation rate ( Barton et al. 2000t Lambas et"ai1 |2003; 



iNikolic et ail 120041: lEllison et al l l2008h and its chemical 



properties dDonzelli & Pastorizal 120001; iMarquez et all 12 002: 
Fabbi ano et al] 12004 IKewlev et al J 120061: iMichel-Dansac et al 
2008: iRupke et all 120081: iKewlev et all 120101: ISol Alonso et al 
2010). Studies have shown that in denser large scale environ- 



ments the star form ation rate of galaxies tends to be lower 
dGomez et a D I2003h . It is often suggested that galaxies in 
clusters are older and have t herefore consum ed the limited gas 
available for star formation dLillv et alJll996t) . or that proximity 
to other g alaxies means that the reservoir of infalling gas must 
be shared dLewis et al. 1 12001 . It is also thought that a dominant 
method of reducing the star formation rate in dense environment 
is through ram pressure stripping that removes the gas envelope 
from field galaxies a s they are accreted to groups and clusters 
(lBaloghetalJl2004bh . 



Proximity to a cluster centre is also known to impact the 
morphology of galaxies and clusters have a greater fraction of 
early type galaxies ; i.e. th e so-called morphology-density rela- 
tion ( Dressier et al. 1 1997b . While the accretion of field galax- 
ies into denser environments may strip gas from the galaxy and 
leave a slowly reddening SO galaxy, it will do little else to al- 
ter the morphology. Morphological transformations are therefore 
attrib uted to gravitational interactions with othe r group mem- 
bers (Moore et al. 1996: IWeinmann et alJ 12006). Denser envi- 
ronments also increase the li kelihood of galaxy mergers. This 
is supported by the findings of lMcGee et alJ d2008l) where an en- 
hancement in the number of asymmetric disks is found in group 
environments (these groups typically have velocity dispersions 
less than 700 km s _1 , smaller than the larger clusters considered 
in many other works). The authors also put forward the intrigu- 
ing conclusion that while the groups exhibit a larger fraction of 
galaxies that are bulge dominated, no evidence is found that the 
group environment has any effect on the bulk properties of the 
disk galaxy population. 



1 



C. G. Few et al.: Properties of simulated Milky Way-mass galaxies 



Mergers have a direct impact on both the star formation his- 
tory of the gal axies and the m etal distri bution. Hydrodynamica l 
simulations bv lHernquistl (1 1989b and lBarnes & Hernquistl fl996) 
show that mergers funnel gas into the central regions of galax- 
ies. This would tend to dilute the gas phase metallicity at 
small radii and trigger centrally concentrated star form ation 
(Ruo ke et alJ2010tlMontuori et al.l2O10UPerez et al.l201 ll) . This 
trend is consistent with observations of flattened metallicity 
gradients in luminous and ultralu minous infrared ga laxies that 
are identified as merge d systems {Rupke et al.ll2008l) and in in- 
teract i ng galaxy pai rs (Elliso n et al. 2008; Michel-Dansac et al. 
2008; Kewley_etalJ|2010|) wherelhe influence on star forma- 
tion rates extends t o projected separations of up to 40 h _I kpc 
( Ellis on et al . 200|. 

Interactions are far more common in denser environments 
and one would expect to see the effect of these interactions im- 
printed on the metallicity of clus ter and group galax ies when 
compared with those in the field. ICooper et alj d2008l) observe 
that not only do members of clusters have greater metallicities 
but also that on average galaxies that are closer to other clus- 
ter members have greater oxygen abundances (an effect of order 
0.05 dex). The conclusion of the work is that metallicity effects 
are not driven by the cluster as a whole but only by the spe- 
cific pro ximity of each galaxy to others; consisten t with past 
findings (Balog h et aT]l2004at iMartuiez etal "" 



The purpose of this work is not to reproduce artificial clones 
of the Local Group but rather to gauge the systematic offset in 
properties between field galaxies and those with a similar de- 
gree of interaction with neighbours such as is encountered be- 
tween the Milky Way and Andromeda. These groups are hence- 
forth termed as "loose groups" to be clear that the systems in 
the loose group sample are not Local Group clones from initial 
conditions designed to reproduce the exact layout of the local 
universe but are chosen from cosmological simulations based 
on certain similarities to the Local Group. In doing so we find 
galaxies with a range of merger histories that nonetheless pro- 
duce disk galaxies. It is hoped that the discrepancy between the 
properties of these two samples will place constraints on future 
simulations and present insight into the failings of field galaxy 
simulations when attempting to recover the properties of Local 
Group galaxies. The method employed in producing this sample 
and the properties of the simulation code are described in $2] In 
§|3] we present an analysis of the galaxy disk fractions, colours, 
magnitudes, metallicity gradients and velocity dispersions and 
discuss these results in §j4] 



2. Method 



The inde- 2.1. Simulations 



pendence of metallicity on large-scale environ ment is perhaps 
refuted by the findings of Ellison et al] d2009f) where a resid- 
ual metallicity-environment effect is found observationally even 
after the dependence on luminosity and colour have been ac- 
coun ted for, however these two resu lts are not entirely irreconcil- 
able. Martinez- Vaa uero et al] (l2009t) select simulated dark mat- 
ter systems based on the mass and circular velocity of haloes, the 
mutual proximity of the haloes and the distance to a halo with 
the same mass as the Virgo cluster. These simulations explore 
the properties of the haloes while relaxing these criteria and find 
that the nearness of massive external haloes is the most signifi- 
cant factor determining the dispersion of Local Group systems: 
the coldness of the Local Group can be attributed mostly to its 
isolation from clusters. 



Recent simulations 
and more advanced 
which ameliorate the traditional 



ulations dRobertson etall 120041: 



have employed higher resolution 
supernova feedback prescriptions 
f ailings of galax y sim- 



Scannapieco et alj 


2009; Sanchez-Blazquez et al.l 




Stinsonetairi2010: 


Rahimi et al. 2010; Brooks et a] 






It is now possible (and prudent) to examine the more subtle 
factors that influence galaxy properties. Given the demonstrable 
difference between galaxies in different (albeit extremely so) 
environments it is reasonable to expect that galaxies in loose 
groups such as our own Local Group may differ from those in 
field. 

At present the majority of literature on this topic focusses 
on constrained simulations based upon cosmological initial con- 
ditions that will purposefully give rise to systems with the 
properties of the Local Group dGottloeber et al.ll20l61 iPeiranil 
l2010tlLibeskind et alj201 UlPeirani et alj2012h or the use of ex- 
treme ly large simulation volumes to search for analogous sys- 
tems fSnaith et al.ll2.Ql ll and references therein). Here we pro- 
vide a complementary approach by identifying Local Group ana- 
logues in a suite of hydrodynamical simulations. In what fol- 
lows, we will explore the hypothesis that simulated isolated field 
galaxies can be considered suitable proxies for Local Group ana- 
logues. 



The galaxies presented here were simulated u sing th e adaptive 
mesh refinement code ramses (v3.01)(Teyssier 2002). RAMSES IS 
a three-dimensional Eulerian hydrodynamical code with an N- 
body particle-mesh scheme to compute self-gravity. The mesh 
automatically refines according to the local particle density in 
addition to a static refinement of nested regions that reduces the 
run-time while maintaining high resolution around the galaxy 
of interest. Detail s of the refinement scheme are described 
by iTevssierl d2002l) . ramses includes density- and metallicity- 
dependent radiative cooling rates, using an i onisation equilib- 
rium with an ultra-violet radiative background (Haardt & Madau 
1996). Gas cells with a density p exceeding a star formation 
threshold of po = 0.1 cm -3 form stars at a rate of p — -p/t*. 
The star formation timescale t+ is itself a function of the den- 
sity and the free-fall ti me through the free param eter, to, as fol- 
lows: t-t, = f (p/po)~ 1/2 dDubois & Tevssierl200 ~8T). We use t = 8 
Gyr, corresponding to a 2% star formation efficiency for po = 
0.1 crrT 3 . We use the kinetic feedback mode of ramses that aims 
to reproduce the energetic and chemical enrichment of SNell 
explosions: After 10 7 years, star particles release some mass, 
momentum and energy into a 2-cell radius feedback-sphere cen- 
tred on the star particle. We characterise the initial mass fraction 
of the star particles with the parameter ?7sn = 10% (which cor- 
responds to the mass fraction of stars contributing supernovae 
feedback) and we do not use any mass loading factor. The en- 
ergy injected into the gas phase is in the form of kinetic energy 
with 100% efficiency (i.e. 10 51 erg SN -1 ). Chemical enrichment 
of the gas is followed through the global metallicity Z using a 
yield of 10%. ramses includes a polytropic equation of state with 
an index of 5/3 to prevent non-physical gas fragmentation in all 
gas cells with an hydrogen density larger than the star formation 
density threshold. The temperature threshold is T± = 10 4 K and 
in the analysis and Table [2] in particular, we will set the temper- 
ature of the high-density gas cells to 10 4 K to account for unre- 
solved cold gas. A detailed descripti on of the star formation an d 
feedback treatments may be found in lDubois & Tevssierl J2008), 
though the v3.01 ramses has a different technical implementa- 
tion of the kinetic feedback. 
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Candidate haloes for this work were identified from dark 
matter simulations and then resimulated with baryonic physics 
and a more refined gri d around regions of in t erest, using 
the same technique as in ISanchez -Blazaue z et al.l (|2009). The 
simulations are conducted in a cosmological framework with 
#0=70 km s -1 Mpc -1 , Q m =0.28, Q A =0.72, Q b =0.045, and 
<t 8 =0.8. Two volumes are used with a size of 20 hr l Mpc and 
24 h~ l Mpc and the maximum refinement achieved (16 levels) 
corresponds to 436 pc and 523 pc respectively. The dark matter 
mass resolution is 5.5xl0 6 M and 9.5xl0 6 M respectively. We 
now describe the selection of our candidate haloes. 



2.2. Environment 

A rich literature exists comparing simulated field disk galax- 
ies with observations o f the Milky Way (e. g. Broo k et al.1 
20051: IS anchez-Blaza uez et al.l 12 009: 



2004; Scannameco et al 



Kobayashi & Nakas atol 1201 ll iHouse et al.l 1201 ll iGuedes et alj 



201 1). The purpose of our work here is to determine whether or 
not environments comparable to those of the Local Group - the 
true environment of the Milky Way - result in any measurable 
characteristics which would call into question this fundamental 
tenet of simulation vs observation comparison. 

The dark matter haloes of the galaxies presented here are 
chosen as follows. Dark matter haloes with virial mass (M v ,,>) in 
the range 5x10' '-10 12 M are considered as loose group candi- 
date haloes. Further selection criteria for the loose group sample 
are such that large groups are excluded; each halo must have 
at least one companion with a comparable virial mass at a dis- 
tance of 500-700 kpc but have no haloes more massive than 
5xl0 12 M within 5 Mpc. Please note that the loose groups are 
not specifically constrained to be pairs, rather they may include 
up to four haloes each with M,,,>>10 n M . 

The field sample consists of galaxies that have no other dark 
matter haloes of mass Myir^xlO 11 M within a 3 Mpc radius. 
These field environment galaxies are far more common than the 
loose groups and those presented here were chosen based upon 
properties such as mass, spin factor and number of mergers such 
that a similar range in each of these can be found in both the 
loose group and field samples. The spin factors of the dark mat- 
ter haloes range from 0.007 to 0.08, this range encompasses the 
majority of dark matter halos and is simply used to exclude out- 
liers. The merger trees of the galaxies are used to select loose 
group and isolated haloes to avoid galaxies (in both loose group 
and field samples) that have any significant mergers after z=l. 
This is done to ensure that disk galaxies are formed and also pre- 
vents the increased likelihood of mergers in denser environments 
from dominating the systematic differences. The full sample of 
galaxies covers a range in virial mass from 10 n -10 12 M and 
includes the field galaxies and all the group galaxies in this mass 
range, i.e. massive satellites in addition to the dominant group 
members. 

There are ten loose group galaxies (taken from three groups) 
and nine field galaxies. The dark matter distributions of a field 
galaxy and a loose group are shown in Figure Q] Each sample is 
roughly divided into two different resolutions corresponding to 
the different cosmological volumes. None of the galaxies have 
passed pericenter with one another as of z=0 and so one should 
not expect dramatic tidal or merger effects. Two of the galax- 
ies do have recent mergers, but for the analysis presented here 
are studied at an earlier timestep (that is unaffected by the z=0 
merger, as described later). The details of each galaxy can be 
found in Table [T] It is well known that the stellar mass fraction 
of simulated galaxies is too high, particularly for older stellar 




Fig. 1. Examples of the different environments considered in this 
work. The dark matter haloes of the field galaxy Selene (top) 
and the loose group galaxies Castor, Pollux, Tyndareus and Zeus 
(bottom) are shown as blue particles (lighter colours correspond 
to denser regions). Images are 4x4 Mpc 2 in size and have a depth 
of 4 Mpc. 



populations. This is true of the galaxies presented here (stellar 
masses are stated in Table [TJ with the stellar-to-total mass frac- 
tion a factor of 2 -3 times too high when counting all mass in 
the virial radius dM andelba um et aH l2006t iMoster et aD 120101: 
iLeaufhaud et al.ll2012l) . We do not believe that this has a dras- 
tic impact on this analysis. Firstly, the issue affects galaxies in- 
dependent of environment, so comparisons of the loose group 
and field samples are not systematically offset by this effect. 
Secondly, the early formation of stars would lead to an overly 
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Fig. 2. Major mergers distributions, the top two rows are loose 
group galaxies and the bottom two rows are field galaxies, a con- 
vention that is followed throughout this work. The magnitude of 
the merger is represented by the virial mass of the host halo di- 
vided by the virial mass of the merging body. As stated in the 
text, mass ratios are calculated at the time the haloes first come 
into contact as the mass associated with merging structures is of- 
ten greatly diminished by the time the disk is disrupted. The ma- 
jority of mergers take place at early times when the galaxy halo 
still has a low mass and assembly is mostly hierarchical. We 
find that only mergers shown above to occur further back than 
4 Gyr have any significant influence on the disk structure and 
that mergers can take several Gyr before they have descended far 
enough into the halo to disrupt the disk. A dashed line shows the 
lower bound of the commonly used definition of major mergers. 



massive spheroid. We will show later that the spheroid mass is 
not extreme, although admittedly the disk stars are too concen- 
trated, i.e. disk scale lengths are too short (see 93.21 ). 

A preliminary study of the metallicity gradients in these 
galaxies has been conducted in comparison with metallicity gra- 
dients fou nd in other simulated ga laxies and in semi-numerical 
models in Pilkingtonet al. (2012) (there these galaxies are re- 
ferred to as Ramses Disk Environment Study or RaDES galax- 
ies); that work demonstrated that when compared with other 
models the RaDES galaxies tend to possess shallower metallicity 
gradients and have a more subtle evolution as a function of time, 
something the authors attribute to the m ore uniform star forma - 
tion profile compared to other models (iPilkington et al.ll2012h . 
This issue is discussed further in 33.21 

2.3. Merger History 

The predominant strength of using cosmological simulations to 
develop this suite of galaxies is that a self-consistent merger tree 
provides the basis of the evolution for each galaxy. This ensures 
the conditions in the environments selected here are consistent 
with the current understanding of cosmology and not biased by 
artificial initial conditions that may affect the results. We begin 



by demonstrating that there is at least a superficial similarity be- 
tween the merger histories of the two samples in spite of the 
environmental differences. 

At each timestep a catalogue of ha loes and subhaloes is cre- 
ated using the Adaptahop algorithm ( Aubert et al.l 12004). The 
halo catalogues are then linked into merger trees for the selected 
haloes using the "m ost massive substructure method" detailed in 
iTweed et al.l (l2009h . Under this formalism, at any given branch 
in the tree the most massive progenitor is considered the parent 
halo. 

It is not however a trivial exercise to define mergers in this 
context, and the following definition of a merger is adopted: a 
merger occurs when an object is identified as a subhalo in a 
given output but not in the previous output. In practice, the sub- 
halo continues to exist as an identifiable structure orbiting the 
host halo for several Gyr. Due to dynamical friction, the subhalo 
gradually sinks into the potential well of the host halo and is 
slowly stripped of mass before dissolving completely, at which 
point all particles are attributed to the host halo. We found that 
in general with the merger time definition above there is a delay 
of up to 4 Gyr between the merger time and any real interaction 
with the disk, i.e., there is a delay time between the dark mat- 
ter merger and the galaxy disk merger. The definition above is 
used because the amount of mass loss before the subhalo dis- 
solves is unpredictable and the mass ratio of subhalo to host is 
a more useful quantity for the purpose of evaluating the mag- 
nitude of the merger. The magnitude and timing of mergers for 
each galaxy are shown in Figure [2] This figure shows that most 
mergers a galaxy experiences occur at early times and that these 
early mergers tend to have mass ratios closer to unity. This is 
partly due to universal expansion reducing the merger likelihood 
and partly to the limited size of haloes at early times limiting the 
mass range and making equal mass mergers more likely. 

The traditional definition of a major-merger, Mh os t/Miub ^ 3, 
has been made significantly more generous (Mh os t/^sub ^ 10) 
in this work due to the low number of 3:1 mergers in the sam- 
ple. The low number of such mergers is due to selection effects, 
these galaxies are chosen from dark matter haloes based partly 
on the merger trees and haloes with many mergers were dis- 
counted as unsuitable for hosting disk galaxies. A commonly 
used metric of the merger history is the time at which the last 
major-merger took place. At this point we note that two of the 
galaxies (identified by the names, Castor and Eos) show obvious 
signs of disturbance at z=0 due to an ongoing merger that makes 
the disk of the galaxy difficult to identify. The presence of merg- 
ers may at first be considered counter to the previously men- 
tioned constraint that the galaxies have no recent major mergers, 
in both cases the merger did not appear in the dark matter-only 
simulation and only became apparent following the inclusion of 
baryons. We remedy this in both cases by analysing the galaxy 
disk at a timestep immediately preceding the disturbance. While 
these galaxies appears in the analysis as a low redshift (z^0.03) 
late-type object, they may also be considered as more local ir- 
regular galaxies if the z=0 output is analysed. We conclude this 
section by stressing that the objective here is not to quantify the 
effect of mergers on disk galaxies but to examine what effect the 
environment might have (i.e. via ambient effects) when the in- 
creased instance of mergers is discounted. As such we point to 
Figure|2]as evidence that the loose group galaxies have the same 
diversity of merger history as do the field galaxies. 
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Table 1. List of the disk galaxies simulated for this work, with their total, dark matter, baryonic, stellar and gaseous masses. The 
spatial resolution of the simulation they originate from and their host environment is also quoted. The properties are calculated at 
z=0 except for Castor and Eos as described in 32.31 



Name 


Environment 


Resolution 


M tot 


M DM 




M stc ll al - 


M g!ls 






(pc) 


(M ) 


(M ) 


(M ) 


(M ) 


(Mo) 


Castor 


loose group 


436. 


1.05xl0 12 


8.70x10" 


1.77x10" 


1.10x10" 


6.76x10'" 


Pollux 


loose group 


436. 


4.23x10" 


3.48x10" 


7.54x10'° 


4.99x10'° 


2.55x10'° 


Zeus 


loose group 


436. 


2.33x10" 


1.97x10" 


3.60x10'° 


2.57x10'° 


1.03x10'° 


Tyndareus 


loose group 


436. 


3.30x10" 


2.82x10" 


4.82x10'° 


3.24x10'° 


1.59x10'° 


Apollo 


loose group 


523. 


8.94x10" 


7.39x10" 


1.55x10" 


1.06x10" 


4.89x10'° 


Artemis 


loose group 


523. 


7.45x10" 


6.46x10" 


9.84x10'° 


5.61x10'° 


4.23x10'° 


Daphne 


loose group 


523. 


3.09x10" 


2.58x10" 


5.14x10'° 


2.90x10'° 


2.24x10'° 


Leto 


loose group 


523. 


2.49x10" 


2.05x10" 


4.38x10'° 


2.92x10'° 


1.46x10'° 


Luke 


loose group 


523. 


1.13X10 12 


9.40x10" 


1.88x10" 


1.07x10" 


8.13x10'° 


Leia 


loose group 


523. 


3.93x10" 


3.25x10" 


6.76x10'° 


4.75x10'° 


2.01x10'° 


Ben 


field 


523. 


7.74x10" 


6.42x10" 


1.32x10" 


8.16x10'" 


5.08x10'" 


Tethys 


field 


523. 


7.21x10" 


5.94x10" 


1.27x10" 


8.21x10'° 


4.51x10'° 


Krios 


field 


523. 


5.68x10" 


4.78x10" 


8.98x10'° 


6.15x10'° 


2.84x10'° 


Atlas 


field 


523. 


6.48x10" 


5.55x10" 


9.35x10'° 


6.15x10'° 


3.20x10'° 


Hyperion 


field 


523. 


1.03xl0 12 


8.58x10" 


1.72x10" 


1.13x10" 


5.92x10'° 


Eos 


field 


436. 


4.64x10" 


3.92x10" 


7.19x10'° 


4.68x10'° 


2.52x10'° 


Helios 


field 


436. 


1.05xl0 12 


8.91x10" 


1.62x10" 


1.23x10" 


3.93x10'° 


Selene 


field 


436. 


6.07x10" 


5.09x10" 


9.87x10'° 


6.71x10'° 


3.16x10'° 


Oceanus 


field 


436. 


1.12xl0 12 


9.19x10" 


1.97x10" 


1.41x10" 


5.59x10'° 



2.4. Disk Decomposition 

We now separate the stellar particles into spheroidal and disk 
popul ations us i ng a kinematic decomposition similar to that 
of lAbadi et alj d2003l) . Details on the decomposition for each 
galaxy can be found in the appendix where the orbital circular- 
ity distributions are presented. In short, stars are assumed to be- 
long to either a spheroidal or a rotating disk component through 
analysis of the orbital circularity, i.e. the ratio of their angular 
momentum, J z to the circular orbit angular momentum 7 c i rc for 
a given particle energy. The successes of this method are shown 
in Figure [AT2l It should be noted that while the J z /J c i TC distri- 
bution has two components, in the intermediate region, particles 
are stochastically attributed to each component. Thus it is not 
necessarily true that stars identified as belonging to the spheroid 
did not form as disk stars. As such the spheroidal component will 
contain disturbed disk stars (which is arguably appropriate) and 
more critically, circularly rotating bulge stars will be attributed 
to the inner regions of the decomposed disk. This does not skew 
the spheroid-to-disk mass ratio as the vast majority of star parti- 
cles have equal mass, furthermore the order of this effect should 
be small as shown in Figure [Aj] which plots the star formation 
history for the entire galaxy (solid line) and the disk (dashed 
line) and demonstrates that the selected disk stars well represent 
the stars formed at late times. 

In the analysis that follows we define a disk annulus to ex- 
clude contamination by bulge stars and avoid halo stars at the 
disk edge. Figure [3] shows the rotation curves for the galaxies 
(and circular velocities for each phase of matter) and choice of 
disk annulus (indicated by diamond symbols at an arbitrary ver- 
tical position). The outer extent of the stellar disk region is first 
constrained by examination of the rotation curve of the young 
stars (less than 100 Myr old), the departure of the young stars 
rotation curve (solid red line) from the gaseous rotation curve 
(solid blue line) is a useful indicator of the stellar disk edge. 
Consideration is also given to the stellar density profiles and 
metallicity gradients of each galaxy, in several cases the den- 
sity profile or metallicity gradient extends beyond or falls short 



of the young disk edge. The final disk annulus is conservative 
to allow for gradients to be measured for each property over 
a consistent region while avoiding bias from unusual features. 
Despite this there are cases where the region over which gradi- 
ents and scale-length are determined has been changed to reflect 
the characteristics of the galaxy in question. A summary of the 
galaxy properties can be found in Tableland star formation his- 
tories for each galaxy are plotted in Figure |A~T1 

2.5. Mock Observations 

The galaxy's stellar and gaseous distributions, ages and metallic- 
ities h ave been used by the ray tracing program Sunrise ( Jonsson 
2006) to produce m ock images. T he Starburst99 stellar popula- 
tion models dLeitherer et alJ ll999) define the colour and magni- 
tude of stellar particles. Scattering and extinction are determined 
by assuming that dust follows the gas phase metallicity distribu- 
tion. Mock images of the galaxies may be found in Figure [4] 
each image being 50x50 kpc 2 in size and produced using SDSS 
g, r and i filters. We draw the reader's attention to the asymme- 
try of the more extended disks of Luke and Oceanus, the warped 
disks of Castor, Tyndareus, Krios and Hyperion and to the red, 
spheroid-dominated Helios. 

3. Results 

3.1. Disk Fraction 

One of the potential differences that may be seen between the 
galaxies in different environments is the bulge-to-disk ratio. 
There is a body of evi dence suggestin g that a morphology- 
density relation exi sts dDresslerj 1980t iGiuricin et alJ |l995; 
Bam ford et al.l |2009|) due to harassment by neighbours and 
the increased likelihood of mergers. The spheroid-to-disk stel- 
lar mass ratios for the sample are shown in Figure [5] and it 
is immediately apparent that the majority are disk dominated 
(spheroid/total<l) in spite of the spheroid containing the net 
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Fig. 3. Rotation curves showing the circular velocity for stars (red), gas (blue) and dark matter (black) as dashed lines and rotation 
velocities (solid lines) of young stars (age<100 Myr at z=0) and gas as a function of radius. Two blue diamonds denote the inner and 
outer "disk radii" at an arbitrary vertical position. These values are chosen from inspection of the rotation curve of the gas, surface 
density maps, stellar surface density profiles and metallicity profiles. The most conservative choice is made to avoid spurious fits 
but to maintain consistency throughout the analysis. 



mass of the halo in addition to the bulge. Only the smallest 
group members are spheroid dominated (Leto, Tyndareus and 
Zeus) with many of the larger galaxies having disk masses ex- 
ceeding spheroid masses by a factor of 2-4. There is no no- 
ticeable tendency for the galaxies with more disturbed disks 
and smaller galaxies to be spheroid dominated; Artemis has a 
disturbed disk and yet still exhibits a spheroid-to-disk ratio of 
0.326, likewise the low mass galaxies Daphne, Pollux and Eos 
have quite strong disk components in contrast with the similarly 
less massive (M v ,,y<4xlO n M ) galaxies Zeus, Tyndareus and 



Leto. Given the typical problem of forming too many stars, par- 
ticularly at early times, in simulated galaxies one might expect 
the spheroid mass to dominate over disk mass but this does not 
appear to be the case; although admittedly the spheroid-to-disk 
mass ratios of t hese galaxies are still not as low as typical bright 
spiral galaxies dWeinzirl et al.l l2009b . Also note from Figure [3] 
that although the characteristic central peak in the rotation curve 
(associated with excessive concentrations of mass) is present, 
the decline of the rotation curve has more to do with the concen- 
tration of the disk. This issue affects all galaxies irrespective of 
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Castor Pollux Zeus Tyndareus Apollo 
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Fig. 4. Mock images of the galaxies. Face- and edge-on views are separated by dashed lines, images are created by combining SDSS 
g, r and i filters and are each 50x50 kpc 2 in size. Simulation outputs at z=0 are used for all galaxies except Castor and Eos where 
the galaxies are shown at z^0.03 
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environment and thus biases in the star formation history do not 
render comparison of the two environment samples unreliable. 

The lack of a clear separation in the loose group and field 
populations in the spheroid-to-disk plot may point to the fact that 
the environments here are not sufficiently different to allow the 
galaxies to manifest different disk properties and that galaxies 
differ significantly only if they inhabit more extreme overdensi- 
ties. It also reflects the dynamics of these groups. None of the 
galaxies have undergone much interaction with any other mas- 
sive group member, having not passed pericenter with one an- 
other at z=0. This removes harassment by massive galaxies as 
a possible source of disruption and leaves only the possibility 
of complete mergers with smaller satellite galaxies than those 
shown in Figure[2]as a possible explanation. 
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Fig. 5. Spheroid and disk masses as determined by kinematic de- 
composition of the z~0 stellar phase. Squares and triangles in- 
dicate the loose group and field populations respectively. Radial 
lines have constant bulge-disk ratios. 



3.2. Metallicity Gradients 

We now examine the metallicity gradients of the galaxies for evi- 
dence of environmental influences. Metallicity grad ients of inter- 
acting and merged galaxies are known to be fl atter (Ellison et al.l 
2008; iRunke et alJlloOSt iKewlev et aUl201(t iPerez et al.ll2oTlh 
and as such, abundance gradients provide a probe of the dy- 
namical mixing. To make comparisons with observed HII gra- 
dients we select stars that are younger than 100 Myr from the 
kinematic decomposition and also employing the spatial con- 
straints described in §2.41 Young stellar metallicity gradients 
are shown in Figure [6] as a function of total galaxy mass. The 
gradients exhibited by the RaDES galaxies range from -0.07 to 
-0.02 dex kpc -1 , consistent with observations by Zaritskv et al.l 
(11994 of spiral galaxies in the field (-0.23 1 to 0.0 21 dex kpc' 1 ) 
Gradie nts are also calculated for spiral galaxies in Ivan Zee et alJ 
(1998) spanning -0.07 to -0.04 dex kpc' 1 and iGarnett et all 



(119971) with a range of -0.083 to -0.020 dex kpc' 1 , the RaDES 
galaxies are remarkably close to these values. 

Metallicity gradients are thought to be flatter for galaxies in 
denser environments and has been demonstrated to be true for 
close interacting binaries bv lKewlev et al.1 ((2010) where HII re- 
gion metallicity gradients are found between -0.040 to -0.007 
dex kpc' 1 . While there is no obvious distinction between the 
two samples presented here it is worth bearing in mind that the 
RaDES loose group galaxies are an order of magni tude more 
distan t from each other than the galaxy pairs in IKewlev et al.l 
d2oToh . Furthermore the simulated gradients are not dramati- 
cally inconsistent with those measured for interacting binaries. 
We note that the simulated loose group galaxies do not have 
appreciably flatter young stellar gradients (with the exception 
of Leid), however they likewise do not have steeper gradients, 
which leaves the possibility that the statistical sizes of the sam- 
ples here may be too small to probe such a slight effect. 
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Fig. 6. Metallicity gradients of the disk stars at z~0 plotted 
against the total mass of the galaxy halo, triangles represent field 
galaxies, squares are loose group galaxies. A trend toward shal- 
lower gradients with increasing mass is evident for both environ- 
ment samples with a slight offset to flatter gradients for the loose 
group galaxies. 

Another feature of Figure [6] is the trend for less massive 
galaxies to have steeper gradients. Leia is a notable outlier from 
this trend, having a particularly flat metallicity distribution, how- 
ever when the gradient is calculated on stars of all ages the trend 
remains and Leia does not appear to be peculiar. This is counter 
to what might be expected since these galaxies have less massive 
dark matter haloes and therefore may be m ore easily perturbed 
and ha ve flattened metallicity gradients. In Prantzos & Boissier 
(2000), cosmologically motivated scaling relations are used to 
demonstrate that the metallicity gradients of spiral galaxies are 
steeper in less massive galaxies when expressed in dex kpc' 1 
but not so when expressed in dex R d 1 where Rd is the disk scale 
length. This is a consequence of the shorter scale-length of the 
less massive disks arising from a steeper star formation rate pro- 
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Table 2. Galaxy properties of the simulated disks at z~0: scale length, metallicity gradient (for stars with age less than 100 Myr), 
cold gas mass-weighted metallicity average (T<1.5xl0 4 K), stellar mass (kinematic ally defined disk stars), cold gas mass, and 
magnitudes (b and r SDSS filters). Note that for the gas mass and metallicity determinations, spatial cuts are used to exclude 
satellites. 
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file which results in a greater metal production rate in the inner 
disk compared with th e disk periphery. This behaviour is sup- 
ported by observations (iGarnett et alJI 19971 Ivan Zee et al . 1998) 
where it is shown that less luminous spiral galaxies have steeper 
gradients than brighter galaxies when the absolute gradient is 
measured. When observed gradients are normalised to the disk 
scale-length however no sign ificant variation wit h luminosity is 
apparent. Another finding of Garn et! et al.l (Q997) is that the dis- 
persion of absolute gradients is larger for less luminous galaxies 
but when normalised the dispersion is consistent with brighter 
galaxies pointing to the existence of some degree of co-evolution 
of metallicity and density profiles. The scale length normalised 
metallicity gradients of young stars are shown in Figure [7] and 
confirm that metallicity gradients likely have a common origin 
wi th the stellar density p rofile. This is explored in more depth 
in Pilking ton et al.l d2012l) where clear links are made between 
the star formation profile and metallicity gradients. The abso- 
lute metallicity gradients of th e RaDES galaxies a re consistent 
with values found in literature (IGarnett et al.ll 1997b however the 
normalised gradients are around an order of magnitude flatter 
suggesting that the measured density profiles are too steep, a 
feature that is consistent with the known issue of excess star 
forma tion at early times and peaked rotation curves in simula- 
tions dNavarro & Ben3ll99lt iGovernato et all 120041: iGuo et alJ 
l2010tlSawala et alj201 lb. Aresult that is relevant to this work is 
that of lDutil&Rovl(ll999l) . Here the authors find that HII metal- 
licity gradients are flatter for earlier morphological types, but 
critically, that the trend is weaker when the gradients are nor- 
malised to some isophotal or effective radius. While we do not 
show the correlation between gradient and morphology (indeed 
no attempt is made to identify the classical morphology of these 
galaxies) it is clear that the metal gradient has some degree of 
co-evolution with the scale length. 
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Fig. 7. Young stellar metallicity gradients at z~0 normalised by 
stellar scale length plotted against the total galaxy mass, symbols 
are as in Figure [6] 



3.3. Colour-Magnitude Diagram 

We use Sunrise to produce seven different projections rang- 
ing from face-on to edge-on and display the values in a 
colour-magnitude diagram (Figure [8) that overplots the ob- 
served (uncorrected) co lour-magnitude diagram from SDSS data 
(Bailin & Harris 2008). We represent the change in magnitude 
and colour as a function of projection angle for each galaxy in 
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Figure[8]with an arc that starts with a symbol denoting the face- 
on projection (symbol type follows paper conventions), the end 
of the arc denotes an edge on projection. Face-on magnitudes 
for each galaxy are given in Table [2] Almost all the galaxies 
populate the blue cloud with only Helios appearing within the 
red sequence. Much as expected, as the galaxies becomes more 
inclined they appear to dim and redden; many of the galaxies 
therefore have edge-on projections that stray into the dimmest 
end of the red sequence. There is a selection effect at work here 
as the galaxies are a priori chosen to be disk galaxies with on- 
going star formation and that only Helios is particularly red is 
reassuring. The colour of the galaxies reflects that while the disk 
stars are too concentrated there is no critical over production of 
stars at high redshift which would bias the galaxies to the red 
regions of Figure [8] 

There is no apparent separation of the two samples once the 
obvious outlier of Helios has been discounted. There is some 
observational evidence that the colour distribution of late-type 
galaxies is only weakly dependent on environment and that it 
is more strongly in fluenced by the lum inosity or mass through 
intrinsic evolution (Bal ogh et al.l l2004b). The distribution in r- 
band magnitude is consistent with the mass of each galaxy and 
the colour of galaxies can be considered a probe of their star for- 
mation history (Figure lATTI ). Taking Helios as an example, we 
see an initially prolonged star formation phase in comparison 
with other galaxies that have star formation rates that are less dis- 
parate throughout time. We examine the impulsive interactions 
of the disk with respect to kinematics in the next subsection. 
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Fig. 8. Colour-magnitude diagram. The background data are 
SDSS galaxies with no inclination correction. Symbols follow 
the paper convention and denote the face-on values. The tails 
traces the change in orientation from face- to edge-on. 



3.4. Sigma-Age Relation 

We now move away from morphology and chemistry to examine 
kinematics and use the temporal behaviour of the stellar velocity 



dispersion (<x) for this purpose. The velocity dispersion of a re- 
gion analogous to the "solar neighbourhood" is used to quantify 
the influence that external interactions have on the kinematics of 
the disk. To remove the bias that arises from the velocity disper- 
sion gradient as a function of disk radius we select stars from 
an annulus of thickness 2 kpc, centred on 3 times the disk scale 
length and a height of less than 3 kpc above and below the equa- 
torial plane. Figure|9]shows the velocity dispersion of stars, cr, as 
a function of their age at z=0. While many of the RaDES galax- 
ies attain the observe d velocity dispersion of the Milky Way (10- 
20 km s" 1 found by ISoubiran & Girardl d2005l) : ISoubiran et alj 
(2008); Holmberg et^alj (|2007j)) when the youngest stars are con- 
sidered, the older populations have far greater dispersions than 
observed. The cr-age dispersion relations shown in Figure [9] ex- 
hibit a greater increase in dispersion as a function of age and the 
appearance of more discrete steps than are apparent in observa- 
tions. Greater velocity dispersions are to be expected in simula- 
tions as less than ten resolution elements are found in the would- 
be thin-disk (if such a structure were resolvable). The high early 
velocity dispersions found even in field galaxies with the fewest 
mergers (e.g. Krios or Selene) may indicate that these galaxies 
experience excess kinematic heating from mergers. Another pos- 
sibility is that very early in the galaxies formation cold flows 
(rather than discrete mergers) are responsible for forming the 
majority of stars. As it is well known that cosmologically sim- 
ulated galaxies form too many stars at early times (see Table Q] 
for the stellar mass fractions) it is possible that the high velocity 
dispersion of old stars is simply reflecting this shortcoming. We 
note however that maps of the gas distribution at high redshift 
clearly show discrete gaseous objects merging and no evidence 
is found for cold streams after z=2. 

The cr-age relations do not at first display any obvious dis- 
tinction between field and loose group environments, so a more 
rigorous analysis is called for. We have attempted to quantify 
the difference between a "stepped" and a smooth cr-age rela- 
tion by looking for spikes in the age-derivative of cr(age). If the 
age-derivative of cr(age) exhibits spikes above some significance 
level (the step-threshold) it will betray the existence of steps in 
the velocity dispersion. Galaxies that have steps with magnitude 
exceeding a threshold (the step-threshold) that is a factor of f3 
greater than the average of the age-derivative of cr are defined 
as having a "stepped profile". The value of /3 is chosen in the 
range 4-7 as values that are too low will not discriminate small 
variations in the age-derivative from larger steps and a value that 
is too large will classify all steps as being normal. This test is 
repeated for different spatial selections in each galaxy (annuli 
with 2-4xR,j in radius with a width of +2-4 kpc including all 
stars within 3 kpc of the disk plane) to avoid bias that may arise 
from sampling a spiral arm and for different values of beta in the 
range 4-7. In each case the field sample has approximately half 
the number of galaxies with step features than does the loose 
group sample. An example is shown in Figure [9] for /3 — 5 (for 
stars in an annulus of thickness 6 kpc centred on 2 x and a 
height of 3 kpc above and below the equatorial plane) in which 
seven out of ten loose group galaxies have stepped profiles while 
only three of the nine field galaxies do. The binning in age usu- 
ally places at least several hundred stars in each age bin making 
the uncertainty in the velocity dispersion (and propagated uncer- 
tainty in the discrete age-derivative) very small, these are plotted 
as 2 sigma error bars in Figure [9] In the extreme case of taking 
/3 = 8 none of the field galaxies have detected steps while four 
of the loose group galaxies do. 

While each choice of annulus and fi results in the loose 
group sample having a greater fraction of galaxies with stepped 
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Fig. 9. Present day (z=0 in most cases, see §2.3\ stellar velocity dispersion as a function of age for an annulus of thickness 6 kpc 
centred on 2 x ^ and a height of 3 kpc above and below the equatorial plane. The lower line is the age-derivative of this function, 
dcr/dage. Horizontal dotted lines define the zero point and a step-threshold that is 5 x (dcr/dage). 



profiles than the field sample, the number of galaxies studied 
here is not enough for this result to be conclusive and a more 
complete sample would be required to make it so. Nevertheless 
we now discuss the implications of the trend if it is real. The 
physical mechanism shaping the cr-age relation can be simply 
thought of as stars being heated to greater dispersion by all 
mergers subsequent to their formatio n and to an extent that de- 
pends on the severity of the merger ( Villalobo s & Helmill2008t 
iDi Matteo et "ai1 l201lh . A series of gentle mergers therefore re- 
sults in a smoother decl ine in dispersion towards younger stars 
dKazantzidis et alj |2008) while a large and disruptive merger ex- 
cites all stars formed previously to a plateau that gives a more 



step-like appearance to the cr-age re lation ([Brook et all 12004). 
This contrasts with the conclusions of lHouse et alj d201 ll) where 
stars are found to form with a velocity dispersion and retain it 
as a signature of the gas state at that time, however in either 
scenario the analysis performed in this work is a valid measure 
of the effect of disk disruption on stellar dynamics. The results 
shown in this work tentatively suggest that in spite of a super- 
ficially similar major-merger history, there may be some dif- 
ferences in the mergers experienced by the galaxy disk or the 
interaction histories. Firstly galaxy disk major mergers could 
impact the galax i es differentially depending on the gas fraction 
(ICoxet al.ll2006t iHopkins et 10120091 lLotz et aLlfeolOl) . the or- 
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bital configuration (Bar nesll2002tlRobertson et alJl2006h and the 
large-scale tidal field (Martig & Bournaudll2008l) . Secondly, the 
number of minor mergers or interactions with orbitin g satellites 
could have an impact in shaping the cr-age rela tion dOuinn et alJ 
l!993HAbadi et alJl2003l:lBournaud et al. 2005). Third, the inter- 
action of the galax y with the intragroup m edium may have an 
impact as found by Bek ki & Couch] (1201 ll) wherein the authors 
find that repetitive harassment in groups can lead to star forma- 
tion bursts and disk heating. 



4. Conclusions 

We present a suite of cosmological simulations with the inten- 
tion of comparing field galaxies with galaxies in Local Group 
environments. The galaxies are taken from cosmological simu- 
lations where a zoom method is used to allow sub-kpc resolution 
while simultaneously accounting for large scale structure forma- 
tion. A kinematic decomposition has been performed to separate 
disk stars from spheroid stars and we have analysed the morphol- 
ogy of the galaxies. While all galaxies studied here have larger 
stellar mass fractions than observations dictate they should, this 
issue is uniform and does not affect comparisons between the 
two environments. We have also examined the metallicity gra- 
dients finding trends with mass but a very weak or non-existent 
correlation with environment. Finally the stellar velocity disper- 
sion is studied and evidence of a dependence on environment is 
found in the signature of impulsive heating in group galaxies. 
The results of this work are pertinent to the comparison of sim- 
ulated field galaxies with observations of the Milky Way. The 
conclusions of this work are summarised here: 

1. No distinction between loose group and field galaxies is seen 
when considering the spheroid-to-disk ratio although exami- 
nation of galaxies with greater spheroidal components shows 
that they all have interactions that disturb their disk rather 
than forming from kinematically hotter gas. However this is 
far from conclusive as there are only four galaxies with sig- 
nificantly higher spheroidal components of the total 19. 

2. Metallicity gradients of loose group galaxies are very simi- 
lar to those of field galaxies with the same disk mass, a result 
that is still c onsistent with obser vations of strongly interact- 
ing galaxies dKewlev et al.ll2010l) though non conclusive ev- 
idence is seen that loose group galaxies should have signif- 
icantly flatter gradients compared with their counterparts in 
the fi eld. The absolute gradients are consistent with observa- 
tions dZaritskv et al.1ll994tlGarnett et al.lll997Hvan Zee et all 
1998) yet when normalised by disk scale-length gradients 
are an order of magnitude flatter than observed suggesting 
that density profiles are too concentrated. Observations also 
show that more massive spiral galaxies have flatter gradients, 
this has previously been matc hed by semi-analytical mod- 
els dPrantzos & Boissierl2 000) using scaling relations but the 
trend has now also been shown to exist for our numerically 
simulated galaxies. We also find a link between metallic- 
ity gradients and stellar density gradients that suggests that 
galaxies in the mass range studied here have similar metal 
gradients when expressed in dex/Rj and that variance in this 
value may be attributable to radial migration or disruptive 
mergers. 

3. Examination of the age-velocity dispersion relation reveals 
that as expected the velocity dispersion of old stars in the 
simulated galaxies is greater than observed for the Milky 
Way disk. Loose group galaxies exhibit more stepped rela- 
tions that suggest mergers/harassment do have a greater im- 



pact on the loose group galaxies than field galaxies. This is 
at odds with the apparent similarities in the major merger 
frequency of the two samples and suggests that the major 
merger history of dark matter haloes may not be an accu- 
rate probe of the galaxy disk merger history. We note how- 
ever that the relatively small sample size means that it is not 
conclusive whether or not this is a real effect depending on 
environment or simply because the loose group galaxies hap- 
pen to have more turbulent formation histories independent 
of environment; certainly the individual formation history of 
each galaxy will impact on its velocity dispersion evolution 
and this would be consistent with the other findings of this 
work. 

The main conclusion to come from this work is that in such 
sparse environments where the galaxies are not directly interact- 
ing galaxies exhibit different properties depending on individual 
merger histories and infall rates but that loose groups environ- 
ments are only very weakly different to the field. It has been 
suggested for cluster galaxies that it is likely that galaxies are 
shaped more by direct mergers and their own secular behaviour 
rather than the large-scale environment that im pacts the afore- 
mentioned only indirectly (McGeeetal. 2008). Structures fur- 
ther than ~ 1 Mpc distant have little influence for galaxies with 
masses presented here. 

The method used to define mergers in this work may be in- 
sufficient to link mergers to the signatures of the impact they 
have on the disk properties and a future study to follow this 
should develop a larger suite using only dark matter simulations 
to quantify the satellite distribution and minor merger rates with 
a greater statistical significance. Simulations at higher resolution 
should also be employed to determine conclusively if any sys- 
tematic difference in metallicity gradients exists. We finish by 
stating that at the resolutions considered here simulated galax- 
ies may be safely compared with Milky Way properties whether 
they inhabit loose group or field environments, however atten- 
tion must be given to the aggregated merger properties, mass 
and internal structure for such comparisons to be meaningful. 
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Appendix A: Additional Galaxy Properties 

Here we describe some properties of the galaxies that are not directly discussed in the main paper body but are still useful as 
supplementary material. We begin with some general comments on particular features of each galaxy that are of interest. 

- Castor is the only galaxy in the sample to exhibit a bar, perhaps reflecting the greater resolution or the more isotropic nature 
of its group (compared with, for example, the filamentary structure of the Apollo group). It also presents the clearest example 
of spiral structure of the entire sample. The spiral structure presents challenges when quantifying the stellar scale length as the 
arms present a bump in the density profile. The young stars present in the arms mean that this is even more pronounced when 
measuring the brightness profile. Castor also has the most pronounced (and abrupt) disk warp, initially this was believed to be 
evidence of poor resolution at the disk edge but investigation has revealed no particularly favoured alignment of the disk warps 
in this sample. Analysis of Castor has been conducted on a snapshot slightly before z=0 due to the irregular shape induced by a 
late merger in this snapshot. 

- Artemis is unusual in that it has a relatively massive dark matter halo (7.45xlO n M ), reasonably low spheroid-to-disk fraction 
(0.32) and flat metallicity profile (-0.0068 dex/kpc), yet its disk scale length is only 1 .87 kpc and is truncated at a radius of 7 kpc. 
There exists a gaseous polar disk (not dense enough to form stars) and yet the vertical velocity dispersion changes very little as 
a function of age. This suggests that the last major merger experienced by Artemis left star-forming gas with a similar velocity 
dispersion to the older stars, an effect not seen in the other galaxies. 

- Leto is the least massive galaxy within the Apollo group and the most spheroid dominated of the galaxies. This spheroid is 
not composed of older stars as in the other galaxies, there is a significant fraction of the spheroid stars that are young. This is 
the result of a low star formation rate at early times compounded by a recent, disruptive event that is evident in the velocity 
dispersion-age relation. 

- Eos undergoes a merger at late times that leaves it with an extremely irregular morphology at the last time step, analysis of this 
galaxy is performed on a snapshot prior to this event. 

- Helios is the most early type galaxy of all. Despite its great mass it is the reddest galaxy and has young stars with around twice 
the vertical velocity dispersion of much of the rest of the sample and a prolonged early star formation episode, this contrasts 
with the lack of an identifiable late merger to result in such a morphology. 

- Selene has few mergers in its history and is one of the most quiescent galaxies, forming the largest disk fraction of all the 
galaxies and presenting definite spiral structure. 

- Oceanus has the greatest stellar mass in the sample (though it is among others with comparable halo masses) and has a rotating 
gaseous disk that extends as far as 40 kpc from the centre. This disk is dense enough to form stars and hence this galaxy has an 
extremely long scale length (6.63 kpc) and one of the flattest metallicity gradients. 

One of the key ways of understanding the formation of galaxies is by examining the star formation history of the different 
components. The distribution of star formation in time tells us a great deal about how the different components of a galaxy 
form. In the course of this work the signature of mergers were found to be identifiable in the bursts of star formation seen in 
Figure IA.1I These bursts can in some cases be associated with steps in the velocity dispersion described in § 13.41 however we 
find that the magnitude of the star formation bursts is a poor indicator of the strength of the kinematic disturbance induced by 
the merger from which they both originate. The dichotomy of the star formation burst and the kinematic effects of mergers is yet 
more evidence that the signature of a merger depends on the gas fraction or phase space configuration of the merging bodies. 
Note the restrained recent star formation of the disrupted galaxy Artemis compared with the more disk dominated Apollo or Oceanus. 

For the analysis of the simulated disk galaxies presented in this work to be as uniform as possible we performed a kinematic 
decomposi tion to examine th e disk stars with minimal contamination from halo and bulge stars. The decomposition employed 
follows the Abadietal ] (120031) method of using the angular momentum of stars compared with the angular momentum expected for 
rotating stars. The distribution of the relative angular momentum is shown in Figure [A2l with blue and red lines highlighting the 
distribution on disk and spheroid stars respectively. 
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Fig. A.l. Star formation histories of the sample galaxies. The star formation rate of all stars within the virial radius at z=0 are shown 
by the solid line, the dashed line is all stars tagged as disk stars at z=0. 
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Fig.A.2. 7 z /y c irc distribution of stars within the virial radius (black). Blue and red lines shows the distribution of the disk and 
spheroid components respectively. Note the existence of a third intermediate component included in the distribution that is associated 
with the disk in some of the galaxies. 
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